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Abstract 



In the context of MSSM, a novel improving procedure based on the renor- 
malization group equation is applied to the effective potential in the Higgs 
sector. We focus on the one- loop radiative corrections computed in Landau 
gauge by using the mass independent renormalization scheme DR. Thanks to 
the decoupling theorem, the well-known multimass scale problem is circum- 
vented by switching to a new effective field theory every time a new particle 
threshold is encountered. We find that, for any field configuration, there is 
a convenient renormalization scale Q* at which the loop expansion respects 
the perturbation series hierarchy and the theory retains the vital property of 
stability. 
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I. INTRODUCTION 



The standard model (SM) forms the bedrock of modern high-energy physics, and ac- 
curately describes all physical phenomena at energy scales up to ~ 100 GeV (electroweak 
scale) . Even in the absence of a grand unification of strong and electroweak (EW) forces at 
a very high energy scale, it is clear that the SM must be modified to incorporate the effects 
of gravity at the Planck scale (Mp ~ 10 19 GeV). In this context, the EW scale is put in 
by hand and there is no idea about its origin: it is a completely free input parameter. In a 
more complete theory one would like to understand its origin in terms of more fundamental 
parameters (e.g., Mp), but such a theory would need more structure than present in SM. 

Moreover, in a general QFT containing an elementary scalar the mass of this scalar would 
be naturally at the scale of the cutoff of the theory^] due to the quadratically divergent loop 
corrections to the scalar self energy (hierarchy problem [p]] ) . If one wants to protect the scalar 
masses from getting these large loop corrections, one needs to introduce a new symmetry. 
The only known such symmetry is supersymmetry (SUSY) || which relates properties of 
bosons and fermions. Although there are various ways in which SUSY might be joined with 
the SM, for simplicity one can pursue a minimal construction and attempt to write down 
a Lagrangian, which is the most general effective Lagrangian invariant under SUSY trans- 
formations up to soft-breaking terms. This then results in the Lagrangian of the Minimal 
Supersymmetric Standard Model (MSSM) ||, where each SM particle is accompanied by a 
super symmetric partner. Supersymmetry also requires two Higgs doublets, as opposed to 
the single Higgs doublet of the SM. 

One of the most remarkable features of the MSSM is that it offers a plausible senario 
for SU(2)l x U(l)y symmetry breaking. However, in this senario one has still to enforce 
phenomenologically a bounded from below potential and the absense of directions in field 
space that may induce a spontaneous breaking of electric and/or color charge symmetries 
0| (a fact that clearly violates experimental observations). 

From the theoretical point of view, in the SM electric and color charge are certainly 
conserved in an automatical way, since the only fundamental scalar field is the Higgs boson, 
a colorless electroweak doublet. The Higgs potential has a continuum of degenerate minima, 
but these are all physically equivalent and one can always define the unbroken U(l)y gen- 
erator to be the electric charge. On the contrary, in SUSY extensions of SM things become 
more complicated. Since the Higgs sector contains at least two Higgs doublets Hi,H 2 , one 
has to check that the minimum of the Higgs potential still occurs for values of Hi, H 2 which 
are appropriately aligned in order to preserve the electric charge. Another perplexity arises 
from the fact that the supersymmetric theory (MSSM) has a large number of additional 
charged and color scalar fields, namely all the sleptons (£) and squarks (q). Consequently, 
conservation of color and electric charge symmetries requires that the minimum of the whole 
potential V(Hi, H 2 , t, q) still occurs at a point in the field space where (q), (£) = (realistic 
or true minimum). 

Yet, the true effective potential in which the vacuum structure is encoded, is a poorly 



1 If the SM were the full story, then the Higgs mass would be naturally of 0(Mp). 
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known object beyond the tree-level approximation. One reason for this is the dependence 
of its loop corrections upon the very many different mass scales present in MSSM, so that 
a renormalization group (RG) analysis becomes rather tricky. In general, when one deals 
with a system possessing a large mass scale Qm, compared with the scale at which one 
discusses physics, large logarithms such as ^(Qm/Qh) always appear which affect the con- 
vergence of the perturbative realization of the potential (loop expansion). In this situation, 
one considers resumming the perturbation series by using the renormalization group equa- 
tion. Nonetheless, in many relistic applications one often has to deal with an additional 
mass scale Q m with the hierarchy <C Q m *C Qm- In MSSM, for instance, one can regard 
Qfi, Qm, Qm as the weak, supersymmetry-breaking and unification scales respectively. When 
we study such a system, we face the problem of multimass scales [f^-|{J: There appear several 
types of logarithms ^(Qm/Qh) and ln(Q m /Q^), while we are able to sum up just a single 
logarithm by means of the RG equation. 

But is it really necessary to take into account these obscure loop corrections? Naively, 
one would argue that they cannot change much of the qualitative pattern of the tree level 
minima. On these grounds, let us recall that classical potential in MSSM receive contribution 
from three sources: D-terms, F-terms and soft-breaking terms. The first of these provides 
the quartic terms Vb = A<£> 4 with A > 0. Now along special directions in field space, known 
as D-flat directions, Vo ~ can occur. If there is a minimum in such a direction, then 
the change in the shape of the potential near the vicinity of that minimum, due to one- 
loop corrections, could create perturbatively a supplementary local minimum lower than 
that already present at tree level.Q In other words, even if the one-loop corrected effective 
potential is expected to differ, point-to-point, only perturbatively from its tree level value, 
it is still possible that its shape be locally modified in such a way that new local minima (or 
at least stationary points) appear. 

Clearly, a carefull RG improvement program is then essential in order to deal with 
the physical problems already mentioned. Trying to illuminate this missing point, in the 



present work a generalized improving procedure based on [10] is applied to the MSSM 



effective potential. The main idea of the method is to make use of the decoupling theorem 



ITJ. By this theorem, it is made sufficient to treat essentially a single log factor at any 
renormalization scale, since all the heavy particles (heavier than that scale) decouple and 
all the light particles (lighter than that scale) yield more or less the same log factors. 

The rest of the paper will be organized as follows. After setting our notation and con- 
ventions, the state of the art concerning the EW symmetry breaking in MSSM is briefly 
reviewed in Sec. II. In the next two sections, we describe the physical difficulties we have 
faced in trying to extrapolate the well known low energy picture of MSSM to higher energies 
and how they have been overcome. A detailed presentation of our numerical implementa- 
tion is given in Sees. V-VI, while Sec. VII explains the reasoning behind our renormalization 
scale choise. The final section is devoted to conclusions and further comments. Detailed 
formulae for the various field dependent MSSM mass eigenstates in a neutral Higgs back- 
ground are presented in Appendix A and some special cases of them in Appendix B. Finally, 



2 Provided that we keep under control the F-term contribution as well. 
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last two Appendices include technical details essential to our numerical procedure, as well 
as a description of the numerical algorithms used for our purposes. 

II. SETTING UP THE FRAME 
A. The Lagrangian 

We are dealing with the MSSM i.e., the simplest supersymmetrized version of the SM. 
The requirements of minimal particle content and matter parity conservation immediately 
dictate the expression of the SU(3)c X SU(2)l x U(1)y invariant superpotential 3 

W = YeVETHieij + Y d Q ja ^fljey + Y u & a U c a H l 2 t l3 + fiH{HUij (2.1) 

where the chiral matter superfields transform as follows 

• Quark Superfields : Q = Q) =: (3, 2, 1/6), U c =: (3, 1, -2/3), D c =: (3, 1, 1/3) 

• Lepton Superfields : L = (j ) =: (1, 2, -1/2), E c =: (1, 1, 1) 

• Higgs Superfields : H x = (f) =: (1, 2, -1/2), tf 2 = (|») =: (1, 2, 1/2) 

under the aforementioned gauge group. Note that the free parameter \i and the 3x3 Yukawa 
matrices Y u , Y^, Y e are generally complex. These ingredients are enough to specify a 
globaly supersymmetric gauge invariant Lagrangian : Csvsy- 

The fact that SUSY is not observed at low energies requires the introduction of ex- 
tra supersymmetry breaking interactions. This breaking, however, should be such that no 
quadratic divergences appear and the technical "solution" to the hierarchy problem is not 
spoiled. Such terms are generally termed "soft" [|TJ] and are interpreted as remnants of the 
spontaneous breaking of local SUSY in the underlying fundamental theory (Supergravity 
(SUGRA) H or Superstrings 0). They include : 

• Mass terms for all scalar fields. 

• Gaugino mass terms. 

• Biliner scalar interactions. 

• Trilinear scalar interactions. 

Consequently, the most general SUSY breaking interaction Lagrangian with real mass terms, 
resulting from spontaneously broken SUGRA in the flat limit (Mp— >oo) is[] 

Csoft = -Imx(BB) - l -M 2 (Wrwr) - \m z (G*GR) 



3 SU(2) indices are denoted by whereas a is a color index and family indices are suppressed. 
Also ei2 = +1. 

r, R are SU(2), 577(3) group indices respectively and again we have supressed family indices. 
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-m^J^I 2 - m\ 2 \H 2 f - m||Q| 2 - m%\Df - m%\V^ 

-m\\L\ 2 - m%\E c \ 2 - (h e D E c + h. A Q ja D c a H\e i j 

+h u Q ja U c a H i 2 e ij + H.c.) - (BfxHiHie i:i + H.c.) (2.2) 

We denote with Hi, H 2 the ordinary Higgs boson doublets and with Q, D c , U c , L, E c the 
squark and slepton scalar fields. On the other hand, the gauginos B, W, G, are considered 
as two component Weyl spinors and h = YA, where A is a 3 x 3 matrix containing the "soft" 
trilinear scalar couplings. All extra soft parameters except masses are generally complex. 
This completes the picture for the low-energy effective theory. 

Altogether one would then need more than 100 real parameters to discribe soft SUSY 
breaking in full generality. Clearly, some simplifying assumptions are necessary if we want 
to achieve something close to a complete study of parameter space. The following set of 
assumptions is adopted: 

• We shall work in the approximation of vanishing intergenerational mixing i.e., 
Y e = diag(r e \ Yl F e 3 ), Y u = diag(Y u 1 , Y 2 , Y u 3 ), Y d = diag(Y d \ Y 2 , Fj 5 ) 
where all non-zero entries are real and positive. 

• fi, B and all trilinear soft couplings are real. 

(The phases of these parameters give large one-loop contributions to CP violating 
quantities, so practically they are quite constrained |0| .) 

• In our analysis we shall also keep Yukawas and trilinear soft couplings from light fam- 
ilies, since their contributions to the one-loop effective potential (our main objective) 
are not always negligible for an arbitrary field configuration. 

A dramatic simplification of the structure of the SUSY-breaking interactions is provided 
either by grand unification assumptions or by superstrings. The simplest possible choise 
arising from coupling the MSSM to minimal N — 1 SUGRA is the following set of assump- 
tions: 

(1) Common gaugino mass : 7711/2 

(2) Common scalar mass : m 

(3) Common trilinear scalar coupling : A 

at a very large scale Mx- More complicated alternatives also exist. However for the time 
being, for the sake of setting our notation, we shall focus on this "universal" senario. 

This reduces the number of free parameters describing SUSY breaking to just four : The 
gaugino mass mi/2, the scalar mass mo, the trilinear and bilinear soft breaking parameters 
A and B. We do assume unification of the gauge couplings at scale Mx — 2 x 10 16 GeV. 
At this point, no specific relation is assumed for the Yukawa couplings at grand unification 
scale (GUT).Q 



5 From a theoretical perspective it would be more natural to impose boundary conditions at the 
(reduced) Planck scale Mp/\^8tt ~ 2.4 x 10 18 GeV or perhaps the string compactification one 
Mq — 5 x 10 17 GeV ]l6|, rather than the GUT scale. In string theory, unification of the gauge 
couplings can be achieved even in the absence of a grand unification. In that case, from our point 
of view, it is of no importance whether we impose the above conditions at Mx or Mq- 
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In order to discuss the physical implications of this Lagrangian at low-energy, we need to 
renormalize the relevant parameters at experimentally accesible energies by computing them 
from a set of coupled renormalization group equations (RGE) [fT7| . (A detailed discussion 
of this procedure will be given in the following sections). This leads us to radiative gauge 
symmetry breaking which we now consider. 



B. RG improvement and electroweak breaking 

The effective potential in general admits a loop expansion |i| as : 



OO 1 

Vef f (KA;Q) = V eff {X a ,0;Q) + £ -0TM(p, = o) (2.3) 



n=l nl 



where V e ff(X a , 0; Q) is a non-trivial function that receives contributions from all orders in 
perturbation theory, X a are the couplings of the theory, Q is a mass parameter upon which 
it explicitly depends and r^ n ^ are the connected proper vertices at zero external momenta. 
Using a particular regularization scheme, one can get rid of the infinities by absorbing 
them into the renormalization of the basic parameters of the theory (A a ,0). The resulting 
parameters become then Q dependent and their evolution is determined by the beta functions 
(3 a (couplings) and the anomalous dimensions 7 (fields). When a similar renormalization 
program is applied to V e ff, order by order in perturbation series, the resulting potential 
function V e // can be cast in a formally Q-dependent loopwise expansion 

00 

V e //(Aa(Q), W); Q) = E V(n) = Vm + V(1) + ■■■ ( 2 - 4 ) 

n=0 

Since Q is an unphysical quantity, the potential in Eq. ( |2.4j) should not depend on the 
choise of it. This property can be achieved, if a change in the scale Q can be compensated 
by appropriate changes in the couplings and field rescalings. Mathematically it can be 
formulated by a RG equation satisfied by the effective potential 

( d d 8 \ 

Q— + - V eff = (2.5) 



dQ d\ a 



which admits the solution 



v e ff = v e ff{\ a {Q),<j>{Q);Q) 

= Q' + Veff : 4> independent shift) 

— + V® + V w + ... (2.6) 

The point is that unless Q' is specifically chosen, then V e ff fails to satisfy the RG equation 
of the usual form 100-0]. The obvious choise for SI' is of course Q' = - £^° =0 V^ n \<p = 0). 
From a practical point of view, adding a field independent piece to V e ff is perfectly harmless 
in problems where only field derivatives of V e ff are of interest and a constant value is assigned 
to Q. 
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An appealing feature of the MSSM is that it can lead to the breaking of electroweak 



symmetry radiatively [17,22]. The correct SU(2)l x U(1)y breaking down to U(l) em is 



achieved by restricting the vacuum expectation values (vevs) on the neutral Higgs manifold 

(Hi)=n 1 5{, (Hi)=H 2 5i (g) = 0, (£} = (2.7) 

Here 5*-, = 1,2) is the well known Kronecker symbol, 7i\, 7i 2 are taken real by gauge 
freedom and the last two equalities have to be satisfied by all the scalar quarks and leptons 
of the model. The low energy classical scalar potential along this direction is then 

2 2 

0°) = mWHtf + m\\H 2 \ 2 + 2ml(H 1 H 2 ) + ^%|^i| 2 - \H 2 \ 2 ) 2 (2.8) 

8 

where m\ = m 2 Hi + /i 2 , m\ = m 2 H2 + /z 2 , m| = fiB and g, g 2 are the U(l) and SU(2) gauge 
couplings.f] 

On the other hand, the one-loop effective Higgs potential of the model, in Landau gauge 
using the DR renormalization scheme [23|j, is 



= k Y, V P W (2-9) 

V 

where k = (167T 2 ) -1 and Ai 2 is the field dependent tree level squared mass matrix of the 
model. We denote by M p the mass eigenvalue of the p th particle and S p , C p are its associated 
spin, color degrees of freedom. M p is the number of its helicity states (p runs over all 
particles), <ft are t ne shifted scalar fields and Q is the renormalization scale. Finally, at 
one-loop order we have 

n' = -VU(\ a (Q),<f>(Q) = 0; Q) (2.10) 



so 



V X - loop = Q' + V (0) + V (1) (2.11) 
The parameters of the potential are taken as running ones, that is they vary with scale 



according to the two loop RGEs with one loop thresholds in all running parameters p4| , |25 

Before carrying on, we would like to comment on a subtle point in definition. We 
know that at any acceptable minimum all the eigenvalues of Ai 2 must be positive, otherwise 
an imaginary part appears in This problem arises because of the non-convexity of 

the tree level potential HPB[, when the gauge symmetry is spontaneously broken. In 



We are using the phase convention \xB < 0, so a TLiH 2 > direction will "deepen" the potential. 
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our case, particularly in the Higgs sector, if we compute the tree level field dependent mass 
eigenvalues for the values (v^v®) which minimize V<® at low -energies, some of them would 
certainly vanish (Goldstone modes). On the other hand, when we try to find the true 
minimum (t>i,t>2) of Vi^i oop , negative mass eigenvalues in the Higgs sector may now appear. 
However, since these negative eigenvalues are of 0(H), we can practically take the absolute 



value of Mp inside the logarithm [28 



Given the low energy scale of EW breaking we must use the renormalization group to 
evolve the parameters of the potential to a convenient scale such as Mz (physical Z-boson 
mass), where the experimental values of the gauge couplings are usually referred. In contrast 
to the tree level potential, Vi_i oop is relatively stable with respect to Q around the electroweak 
scale [p0jp8H30t[. Therefore, the exact scale at which to minimize is no longer critical as long 



as it is in the electroweak range. 

A glance at the superpotential Q2.1|) shows that H 2 couples to the t (s) quarks, while 
Hi only couples to b (s) quarks and r (s)leptons by means of the Yukawa couplings. These 
Yukawa couplings enter the RGEs for square scalar masses with positive sign. So their effect 
is always to decrease the Higgs masses as one evolves the RG equations downward from 
GUT scale to Mz- Combined with the experimental fact that top quark is heavy, this can 
cause the RG-evolved rn 2 H2 to run negative near the EW scale, helping to destabilize the 
point Hi = H2 = and so triggering the breaking of gauge symmetry. If we define 



m i = m 1 + = 1) 2) where 

* 2% art. 

k ^ , ^, „ v „ 1 <9M 2 



| £ (-l) 2S ^2S p + l)C p Af p - f(M 2 p ;Q) (2.12) 



with 




f{x,Q)=x[]n^-l) (2.13) 



then the minimization of the potential yields the following two conditions^ among its pa- 
rameters (all parameters are Q dependent) 

sin 2/3 = - _^ Bfi _, (2.14a) 
1 9 / , , x — ml tan 2 8 

2 m *< M *> = (2 ' 14b) 

where m\{Mz) is the running mass of the Z boson and tan 8 = v^jv\. Note that in the above 
relations include contributions from all particles (see Appendix A), Q takes a constant 



7 For an analytic study of these conditions in the Higgs sector see Ref. [27]. 
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value (Mz) and derivatives are taken with respect to running fields Hi(Q), H.2(Q), so there 
is no contribution from vacuum energy (dQ'/&Hi(Q) = 0). Moreover, since all the scalar 
fields (0) are multiplicative renormalizable, we conclude that if <f){Qo) = 0, for some Qq, 
then 4>{Q) — for every Q. This is what actually happens for all scalar fields except Higgses 
(seeEq. O). 

Of course, we do not merely want SU(2)l x U(l)y to be broken; we also want the Z 
boson to have the experimentally measured mass (Mz)- Furthermore, in view of the strong 
dependence of some weak-scale quantities on tan/?, it is often more convenient to treat Mz 
as an independent input parameter. Since m| = | (g 2 + #1) (^l + fixing mz (equivalently 
Mz) and tan/3 determines both vevs Vi,v 2 - 

For all that, destabilizing the origin by a negative (mass) 2 parameter is not enough to 
ensure the viability of the MSSM one loop scalar potential. We must also make sure that 
the potential is bounded from below for arbitrarily large values of the scalar fields, so that 
Eq. (|2.11 ) will really have a minimum. 



III. ATTEMPTS TO DEAL WITH HIGH ENERGIES 

Extending this well known low-energy picture to high-energies (large field values), one 
is confronted with peculiar effective potential configurations. The simplest generalization is 
allowing Higgs fields, in Eqs. fl2.8|) -( ^9|) , to take arbitrary values keeping at the same time 
the renormalization scale fixed at Mz- However, this assumption leads to an unbounded from 
below (UFB) potential. This realization is clearly physically undesirable. Before explaining 
the reason behind such a failure, let us see what causes this fake instability. 

Let 7ix — Xi, H2 = x 2 and calculate all the field dependent mass eigenvalues in polar 
coordinates 

Xi=rcos9, x 2 = rsm9 (3.1) 

With no harm of generality, as shortly will be seen, we'll only take contributions from 3 rd 
family. We intent to write V\-i oop for r ^> 1. For our approximation to hold, coefficients 
that multiply powers of r should not be arbitrarily small. So, for each r, our approximation 
is valid only for those values of 9 that respect the above constraint. Choosing 9 = n/2 and 
using the formulae in Appendix B the potential becomes after some algebra 



where 



V^ioop = V {0) + V {1) + (field independent piece) (3.2a) 
V i0 \r » 1) = my + l±A r ± ~ t±A T ± (3.2b) 



ii 



64vrV«(r > 1) = r 4 \As(r) +J2 d i U ! ln W ) ( 3 - 2c ) 



8=1 



A = ^ + ^ 1)2 - 3F,V + 9l) + E dtf (3.2d) 



i=i 
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In the above, s(r) = \n(Y 2 r 2 / Q 2 ), Q = Qe 3 / 4 and d h Ui are given in Appendix B. The 
crucial term for determining the behavior of the function presented above is the logarithmic 
coefficient (A). This number due to large top Yukawa coupling is negative, so the whole 
function is UFB. (Note that should we have taken contributions from light families, nothing 
would have changed since top quark Yukawa coupling still dominates). 

Apparently, the main tool of our discussion is the loop expansion. So ultimately one has 
to justify the convergence of loop expansion at high energies, ensuring in this way that only 
the first terms in the series should be kept. Eventually, one has to study the logarithmic 
structure of the effective potential. It has been shown in J7,19] that the L-loop effective 



potential contains logarithmic factors, only up to L-th power, whose magnitude control the 
convergence of the loop expansion. These factors have the general form 

s = A In q 2 (3.3) 

where A is some coupling of the theory, M is a field dependent mass eigenvalue in the 
presence of the background fileds <p c and Q is the renormalization scale. In our case (MSSM), 
large field values generate large (field dependent) mass eigenvalues and obviously large log- 
factors (Q is fixed at M z ). So higher order corrections (2-loops etc.) become significant and 
should also be taken into account. These corrections should raise the potential because the 
theory should be stable and can not have an UFB potential. In conclusion, the one-loop 
approximation to the effective potential renormalized at Mz for large field values is not 
reliable and a different scale choise to control large logs is needed. 

An alternative way of thinking about RG improvement is to view it as a reorganization 
of the perturbation series, in which the first term is the sum of all the leading logarithms, the 
second term represents the next-to-leading logarithms, and so on . The leading log- 

arithms are terms of the form ln(Mi(0)/Q) ln(M2(0) / 'Q) ■ • Tn(M n (0)/Q) and represent 
the most "dangerous" logarithmic terms at each order in perturbation theory.| If we could 
fully sum the multiscale leading and next-to-leading logs, we should have an approximation 
that is useful dispite the existence of widely differing scales. However, when there is more 
than one mass scales present it is less clear how to proceed; no choise of Q will destroy all 
the logarithms simoultaneously. 

This point of view is adopted in the improvement prescription presented in J/]. Let us 
briefly describe it. Imagine a theory (like MSSM) having many coupling constants and mass 
scales. The problem in such multiscale cases is which log-factor we should choose to sum 
up the leading, next-to-leading etc. logs. The best choise would be to take a particle whose 
coupling constant is the largest .0 Calling that particle by label fc = 0we have 



8 To compute the coefficient we need graphs with n-loops. Furthermore, the next-to-leading 
logarithms are proportional to u^ n+l ^ ln(Mi(<^>) / Q) ln(M2(0)/Q) • • • \n(M n {4>) / Q) and the tree level 
potential is counted as a leading logarithm. 

9 In MSSM this role is played by the top quark. 
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s 



AoM M o 



Q' 2 



(3.4a) 



M|(0 C ) X k M fc 2 (0 e ) 
s fc = A fc In — — — = —s + u k where u k = X k In (3.4b) 
Q 2 A M O 2 (0 C ) 

If it happens all the masses M|(0 C ) be of the same order independently of the background 
scalar fields C , then u k variables remain always of order 0(X k ) ^ 0(Ao). So all the log 
factors essentially can be treated as s or (Afc/A )s , since the differences u k are of order 
higher than the first Sq term. In this case, the field dependent scale Q 2 = Mq (0 C ) <^ sq = 
correctly sums all logs in the leading-log series expansion, improving thus succesfully the 
effective potential. 

We have numerically applied the prescription described above to MSSM and suprisingly 
found that with this scale prescription the effective potential values near x\ = ±X2 direction 
(H-diag) are quite unexpected. H-diag direction, for large field values, develops saddle 
points (near 200 TeV! as shown in Fig. la) which lead to an unphysical (UFB) potential. 
Unfortunately, near H-diag direction field dependent mass eigenvalues far away from origin 
are not all of the same order (see Appendix B). So some of the u k are large and one has to 
keep any higher powers of them in the leading-log series expansion. This in turn implies that 
higher loop corrections are significant and can not be neglected. In other words, the absense 
of a unique scale choise (even field dependent), that eliminates large logs to all orders, makes 
the improvement prescription of Ref. for MSSM unreliable. 

This deceptive deadlock of course steams from our careless treatment of different mass 
scales by a single scale parameter Q. Hiding all the heavy particle loop contributions in the 
redefinition of the low energy theory parameters, we can still solve the same RG equation 
for the effective potential by using different effective field theories. This will be the subject 
of the next section. 



IV. AN ALTERNATIVE SENARIO : V-THRESHOLDS 

Recently, the authors of Ref. [|1(| have proposed a nice method to realize the attractive 
conjecture just mentioned. Specifically, one should use the decoupling theorem to handle 
with the problem of many scales in the effective potential, using as decoupling scale for the 
mass eigenvalue M 2 (4>) the scale 

Q 2 = |M 2 (0)| (4.1) 

(recall that Q = Qe 3 ^ 4 ). In other words, we use the following expansion as our impoved 
potenial 

V W = k J2vl 1] 6 l (4.2a) 



0. = 0(Q2_| M 2 ((/)) |) 
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(4.2b) 



= i -^f 1 (2S t + 1) dMi M${$) In (4.2c) 

Despite the appearance of Heaviside functions the effective potential is a continuous function 
(a discontinous potential is physically meaningless). A discontinuity is likely to appear only 
at a decoupling point (threshold). Let's examine what is happening when Q 2 approaches 
the £ th threshold. We have 

• Just above the threshold : 

i=i 

• Just below the threshold : 

i=l 

• At the threshold : 

[V% =0 = fc£ [vf^] £=0 + k [V e W 9 £ ] e=0 
1=1 



But V e {1) oc Mf{4>) ln( | M|(0) | /Q 2 ) and as Q 2 -> |M|(0)| (i.e., e -> 0), then [V^ (1) ^] e=0 = 
and [^ 1 ' ) ^] (E+ = so, the potential function is continuous as it should be. 
On the other hand, defining 

£I = ^]Ct-^]pm (4-3) 

the corrections Ej of Eq. (2.12) become in this new approach 

= E; -^E^HO 2 - |M fe 2 (0)|) ^1 (4.4) 

The last term is obviously zero, so the stationary conditions have the well known form of 
Eqs. (2.14a)-(2.14b) replacing E, with the new E* . 

Let us recall that we are working in the one-loop approximation, which is scale inde- 
pendent up to 2-loop corrections. These corrections, however, may be very large due to the 
presence of potentially large logarithms in V^ n>1 \ Thus, one has to be careful in choosing a 
value of Q which gives the finest approximation to Vi_; oop . We'll come back to this critical 
issue during the numerical method's description, but for completeness we stress here that in 
the above stated framework the potential is indeed bounded from below as Fig. lb shows. 
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V. NUMERICAL PROCEDURE FOR RADIATIVE EW BREAKING 



In this section, following the analysis of [24j] we shall present the procedure we use to 



restrict the SUSY parameter space to those regions, where the requirement of radiative 
breaking of electroweak sector can be correctly implemented. 

The generic problem, we have to tackle here, could be formulated as a set of coupled 
first-order ordinary differential equations (ODEs) for the dependent variables (running pa- 
rameters) t/i, i = 1, 2, . . . , n having the form 

y'M = Mvi,V2,---,y n ) (5.1) 

where the functions /j on the right-hand side are known and t = 21n(Mx/Q) is the inde- 
pendent variable (scale). However, a problem involving ODEs is not completely specified by 
its equations. Initial conditions (values of the running parameters y,) must also be supplied 
at some starting point t a which in our case is: t a = i.e., Q = Mx- The underlying idea 
of any numerical routine, that solves an initial value problem, is to propagate the solution 
Hi over an interval [t a ,t b ] with t b ^ t a , using a step-by-step calculation, which approximates 
values of the solution at discrete specified points in that interval. As far as the MSSM case 
is concerned, we have used special routines from NAG fortran library (which is available to 
us) and particularly an Adams variable-order variale-step algorithm, which provides a choise 



of automatic error control and an option of a sophisticated root-finding technique |3T| . 

The initial conditions at Mx scale will be choosen as simple as possible postponing for 
elsewhere the study of more complicated alternatives. Thus, at the unification point taken 
to be O(10 16 ) GeV, we shall take 

rriQ(Mx) = m bc (M x ) = m 0c (M x ) = m (5.2a) 

m L (M x ) = mzc(Mx) = m (5.2b) 

m Hl (M x ) = m H2 (M x ) = m (5.2c) 

and 

M X {M X ) = M 2 (M X ) = M 3 (M X ) = m 1/2 (5.2d) 

In addition we take equal cubic couplings at M x , i.e., 

A e (M x ) = A d (M x ) = A U (M X ) = A (5.2e) 

Obviously, all our "universal" boundary conditions are family blind. 

Contact has to be made between the low-energy and high energy parameters of the 
theory P^fl. This will be achieved by integrating the RGEs from a superheavy scale, taken 
to be in the neighborhood of 10 16 GeV, down to a scale Qo ^ Mz- In order to make a run 
starting from Mx we must have the values of all parameters at this scale, something difficult 
to achieve. The problem is that the values of many parameters are known experimentally 
at low scales. Also 1-loop stationary contitions are applied at such low energies. However, 
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the values of other parameters, such as soft breaking terms, are most easily understood at 
higher energies, where theoretical simplification (e.g., universality) may be invoked. Thus, 
there is no scale at which both theoretical simplicity and experimental confirmation coexist. 

Moreover, MSSM is characterized by two kinds of parameters, constrained and free. The 
former are constrained by experiment (gauge couplings, quark masses etc.) or by relations 
among themselves (stationary conditions at Mz)- The latter, given the present experimental 
data, can not be constrained by these two critiria and must be viewed as input parameters. 
These are A , B , fi (values at M x ), ^1/2, mo, tan (3 (ratio of vevs at M z ) and M t (physical 
mass for the top quark). In the present framework, B and /io will be determined using the 
numerical integration routines in conjuction with the minimization of the one loop effective 
potential at Mz (i.e., Eqs. (2.14a)-(2.14b)). Namely, stationary conditions^ at Mz will 
give B(Mz), [i{Mz) and then we run back up to Mx in order to determine B and /io- Note 
that the sign of /i is not determined from (2.14a)-(2.14b), thus we must make a choise for 
it. So ultimately the free parameters of the theory are taken to be 

A , m , mi/2, tan/3, M t , sign(/x) (5.3) 

On the other hand, the exact values of the constrained parameters at Mx are affected by the 
choise of values for the free parameters at that scale, since the evolution of all parameters 
are coupled. 

We begin our numerical procedure at the EW scale which we take to be Mz- This is 
an obvious choise, since many experimental quantities are now available at that scale. At 
Q = Mz we take as input the well measured values of the Z mass Mz = 91.1867 ± 0.0020 
GeV, the electromagnetic coupling constant oi'^Mz^jjg = 127.88 ±0.09 |33| and the weak 
mixing angle (sin 2 MIms = 0-2316 - 0.88 x 10~ 7 (M 2 - 160 2 )GeV~ 2 0. In addition to the 
values of the gauge couplings^] at Q = M z one also needs the Yukawa couplings of quarks and 
leptons there. In order to determine (Y T , Y b ) = (V e 3 , V 6 3 ) at M z , we run the gauge couplings 
a and 0J3 from their experimental values^ at Mz down to the b- and r-mass scalesQ using 
three-loop QCD and two- loop QED RGEs ||35|| . The relevant procedure is described with all 
details in |24|| . Finally, the top Yukawa coupling (Y t ) is obtained from the running t-quark 
mass, which for our input value M t = 174 GeV is m t (M z ) ^ M t / (l + ^) ~ 162 GeV. 



10 We emphasize that S^, E?; involved in (2.14) contain contributions from all mass eigenvalues, 
even light ones. These mass eigenvalues are obtained using everywhere analytic formulae except 
for the neutralino sector where numerical diagonalization is performed (Appendix A). 

n The modified minimal subtraction (MS) values for the gauge couplings are related to their 
dimensional reduction (DR) ones through the relations g-^ = — Igfr) where C = 0,2,3 

respectively for the three gauge groups. 

12 We start with the MS values for the gauge couplings at Mz giving a trial input value for the 
strong coupling constant 03 in the vicinity of 0.120. 

13 For the b quark and r lepton our inputs are their physical masses: M b = 5.0 GeV and 
M T = 1.777 GeV. 
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The evolution of all couplings from Mz running upwards to higher energies determines 
the unification scale Mx and the value of the unification coupling ax (in DR) by 

di(Mx) = a 2 (M x ) = a x (5.4) 

Running down from Mx to Mz the trial input value for a 3 has now changed. The entire 
procedure is repeated several times. After one third of iterations is completed, we introduce 
finite part contributions to our low energy inputs such as QCD corrections to the top quark, 
gluino mass and weak mixing angle sm 2 OwIdr- Also the gauge couplings values &i, a 2 (in 
DR) are now determined in terms of Fermi Constant Gf = 1.16639 x 10~ 5 GeV -2 , the Z 
boson mass and the electromagnetic coupling = 137.036 |33]. For details on these 



matters see Ref. [35]. The running t, b, r quark masses at Mz are now defined through the 

m s [Mz) =M f + U f {Q = M z ) (5.5) 

where Mf is the physical mass and Uf is the self energy corrections of the fermion. This new 
procedure is now iterated until convergence is reached. When convergence is established, 
the following actions are taken place at renormalization point Mz : 



• Light Yukawa couplings are being computed dividing the fermion masses |J3| by the 
appropriate Higgs vev. These couplings are considered as fixed with respect to QJ^\ 
On the other hand, the light trilinear scalar couplings do run. 

• The whole set of the MSSM running parameters values (gauge - Yukawa couplings, 
soft masses - couplings and Higgs vevs) is stored for later use (improvement of effective 
potential around EW scale). 

• MSSM running parameters are evolved from their current values at Mx to a lower 
scale Qhigh- Qhigh is taken bigger than the largest field dependent mass eigenvalue 
appearing in V^ l \ when the background scalar fields are restricted in a hypercube of 
side Soo. This point will be clarified later, but for the moment we simply state that 
Qhigh — 55oo is a more or less satisfactory approximation for our purposes. 

Finally, our results should be reconciled with 

• Basic experimental constraints on supersymmetric masses as well as Higgs bosons 
masses as shown in Table I fl33| . 

• Radiative electroweak breaking (existence of real non-zero vevs) . 

• Positivity of eigenvalues of the mass matrix. Negative eigenvalues mean that expanding 
around realistic minimum tachyonic states appear making the true minimum of the 
theory unstable. 

VI. THRESHOLD PARAMETRIZATION FOR /3-FUNCTIONS 

In the minimall low energy SUGRA model being considered, the super-particle spectrum 
is no longer degenerate. This should lead to various course corrections, each one occuring at 



14 Since Yught ~ then = (• • -)Y (Ag ht - thus Y iight {t) ~ const. 
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the super-particle mass thresholds. So the renormalization group /3-functions must be cast 
in a new form, which makes the implementation of the threshold effects evident. Since the 
DR RGEs are mass independent, each super-particle mass determines a boundary between 
two effective theories. Above a particular mass threshold the associated particle is present, 
whereas in the effective theory below the threshold the particle is absent. 

The simplest way to incorporate this is to treat the thresholds as steps in the particle 
content of the renormalization group /3-functions JMRDf. Let's briefly outline the proce- 



dure. Assume that b is the beta function of a running parameter in the DR scheme. The 
corresponding RGE should be integrated from a superlarge scale Mx down to any desirable 
value of Q. As we come down from Mx, as long as we are at scales larger than the heaviest 
particle in the spectrum, we include in b contributions from all particles in the MSSM. When 
we cross the heaviest particle threshold, we switch in a new effective field theory with the 
heaviest particle integrated out and of course a new b. Coming further down in energy, we 
encounter the next particle threshold at which point we switch again to a new effective field 
theory with the two heaviest particles integrated out and a new b. That procedure goes on 
until all particles are exchausted. 

Crossing a particle's threshold means that the renormalization scale has become smaller 
than its physical mass. Hence, we need a condition to determine the exact point of decoupling 
(i.e., decoupling scale). For field configurations in the low energy regime (< 300 GeV) this 
is simply 

Q 2 = m 2 (Q), Q = Qe 3/4 (6.1) 

where m 2 (Q) is the running soft parameter corresponding to the particle.^ Consequently, 
the step functions in RGEs will have the form 

e m = e(Q 2 -m 2 (g)) (6.2) 

Alternatively, for all other field configurations we shall use a different condition to fix the 
decoupling point of a particle. Our condition now is 

Q 2 = \M 2 (<j ) ;Q)\ (6.3) 

where M 2 is the field dependent mass eigenvalue for that particle. Analogously, the step 
functions in RGEs will become 

#Af = #(Q 2 -|M 2 (0;Q)|) (6.4) 



This procedure is generally more accurate than the approximation stated in Eq. ( |6.1|) , but 
in the case of true minimum it introduces non-trivial field dependence, through mass eigen- 
values, in £*, £2 an d the simple stationary conditions (2.14) become extremely involved. 
(For a presentation of the relevant m 2 , M 2 see Appendix C.) 



15 We use the factor e 3 / 4 for compatibility with the "analogous" decoupling in effective potential. 



16 



VII. CHOOSING THE SCALE 



Our starting point is that the full effective potential is independent of the renormalization 
point Q and thus satisfies the RGE 



dV f 



eff 



dQ 







(7.1) 



One can imediately write down the general solution to (7.1) introducing the "running" 
distance t from the initial values as : 



Veff = V eff {\ a {t), <f>{t); Q(t)) = n' + V® + V« + Q(h 2 



(7.2a) 



Q(t) = M x e 



-t/2 



-^ (1) (0 = O) 



(7.2b) 



where \ a (t) are all dimensionless and dimensionfull couplings of the MSSM and <p(t) = ((t)<j) c 
are the running; fields with ((t) = exp{- /* dt'^t')}, 7 (t) beine; the anomalous dimension of 
the field. The key to the usefulness of the RG is that we can choose a value of t such that 
the perturbation series converges more rapidly than the series for t = 0. Moreover, there is 
nothing to stop us choosing a different value of t for each value of <j). 

In order to validate the use of 1-loop effective potential one must ensure that not only the 
couplings be perturbative, but that the loop expansion be convergent as well. In problems 
with only one mass scale RG improvement is straightforward. But for the cases of interest 
here there are several mass scales, so one must think up an improvement prescription. More- 
over, the lack of analytic formulae describing the scale dependence of the quantities involved, 
as well as the absence of any profound physical reasoning for choosing the appropriate scale 
make this task quite obscure. 



Earlier attempts |20| , |28| - p0| , p7| can not offer a substantial aid, since their object was an 
improvement in the low energy region. For example, Ref. PSJ argues that there is a scale Q 



where one-loop stationary configuration coincides with the tree level one. Thus by definition 



(7.3) 



The above definition is non-trivial to implement, so one should approximate Q with an 
average of dominant field dependent masses. For \Hj\ < 10 5 GeV this is a legitimate ap- 
proximation, but when extending for \Hi\ > 10 5 GeV the potential develops an unphysical 
UFB escape along H-diag direction which forces us to search for something else. 

Recently, another point of view have been introduced in Ref. ||10|| . Namely, one should 
compute at a scale Qr where both 



Qr 







(7.4a) 



dVi. 



loop 



dQ 



(7.4b) 



17 



are satisfied. Clearly, for this Qr the radiative corrections to tree level are small and our 
approximation to the full potential has the least Q-dependence. We have tried to implement 
the above prescription to MSSM, but the result is rather discouraging. The field dependent 
spectrum of MSSM always consists of masses at the order of Mz (coming from the neutralino 
sector). So application of ([7.4a]) will lead to a critical scale around Mz, which as we saw 
gives a problematic potential.^ On the other hand, ([7.4b|) seems much more reasonable 
in the sense that the potential will be more scale independent. However, in this case our 
numerical routines can not provide us with a scale for every field point (Hi,H2)- Due to 
the step functions in RGEs the integration routine has to take too many intrinsic steps to 
preserve stability and accuracy of the solution. This in turn affects the performance and 
the numerical root finding facility when invoked for the equation (|7.4b|) . 

To overcome this ambiguity in MSSM, we shall make a physically motivated choise of 
scale which reproduces a well behaved bounded from below effective potential one expects in 
a stable theory. Before carrying on, for clarity reasons, let us sketch briefly the behaviour of 
near H-diag direction. Without harm of generality, we choose to examine the dominant 
contribution of top-stop sector. The corresponding partial sums of are 

\M 2 \ \M 2 I 

Top : P t = -2 CMf In L^ii 9 t , Stops : P T . = CM*. In 6 T . (7.5) 

Q 2 Q 2 

where C is some constant factor and 9 t , 9 T . (i = 1, 2) are the associated step functions. For 
non-zero partial sums we must take Q > \M t \, |MyJ, i.e., In < for all M involved 
(Pt > and Pt % < 0). As we approach H-diag bottom points, the field dependent top-stop 
eigenvalues conspire to produce an extremely largeQ negative contribution to which is 
responsible for the UFB escape. 

However, such a picture can not be reconciled with a perturbation series hierarchy, neither 
with the notion of stability every acceptable physical theory should have. Consequently, Q 
should be taken such as top-stop and similar "heavy" pairs are decoupled, see Figs. 2a-b, 
making loop expansion meaningfull (|^^°- ) | ^ l^ 1 ^)- One such rather conservative choise 

is Q* = 1(T 3 sjn\ + H\ for large Hi (Hi > 10 8 GeV).0 Note that for this choise as H U H 2 
approach Mz, Q* will become lower than Mz making the RGE evolution ambiguous. The 
situation can be improved by defining 

Q* = lO^y?*? + Hl (7.6) 



in 



Of course this is not a real problem, since one can interpret (7.4a) as an order of magnitude 



relation between zero and first orders in loop expansion. So it is better to require -C \ V^\ 

at Qr. However, now there are several scale choises satisfying the previous condition. 

17 Compared to y (0) . 

18 Dependence on Hi, Hz is for practical reasons. Q* should be a smoothly continous function for 
the various field configurations and the only available "free" variables are these. 
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1 f-f^ I 7-^2 

where x = £ log 2 is the order of magnitude of a generalized "radious" in the field space 



ij+v 2 

2 L 2 

and L = 1 GeV makes the log argument dimensionless. 

In the context of the approach we have followed, it seems rather non trivial to rigorously 
define uj(x). So one has to look for a qualitative fixing. Specifically, as x decreases, field 
dependent masses decrease too, therefore Q* should decrease, otherwise the non-zero partial 
sums in could produce points deeper than the EW minimum. Using the cubic spline 
interpolation method H^J^SJ we can give an ansatz for to{x), see Figs. 3a-b, compatible with 
continuity and the following constraints : (a) Q* > M Zl (b) |V^| > \V^\, (c) u(9) ~ —3 in 
order to recover the previously stated conservative choise for Q*, (d) as Hi — > Mz, Q* — > Mz 
i.e., 



UJ[X) 



' UJ C — X 
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where by definition x c 
below 



.4 



log(300^), 



\ 



-0.691121 -4.17735 -4. 

0.389965 -6.31201 -4.2182 
-0.176436 -4.46928 37.1171 

0.115471 -10.7737 27.1562 
-0.0417777 -6.09931 183.38 



(7.7) 



and the coefficients A are shown 

-2.39498 \ 
-4. 

-11.5691 (7.8) 

-10.2237 

-26.1999 j 



To complete this picture one also needs some "boundary scale" Qhigh, which shall provide 
the starting values of the running parameters for the evolution at Q*. One convenient choise 
is an intermediate scale higher than the largest field dependent mass eigenvalue at the 
current field point. Specifically, at some field point TC = max(7i 1 ,7i 2 ) we approximately 
have max fc { J \ Mf (7-^, H 2 )\} ~ H. Let denotes the upper bound order of magnitude of 
the allowed values for the fields. Then valid choises for Qhiqh are : Qhigh ~ Soo =>■ Qf ■ > 



'high 



Soo^l 4 — 2.12 Soo. This is reasonable since the SM is supposed to be a low-energy effective 
theory coming from a more fundamental theory. 

Specifying the scale Qhigh is not enough. We also need to know the values of the running 
parameters and fields there. Since Qhigh is above all thresholds, the required values sould 
not depend on the background fields and the reasonable choise at an arbitrary field point is 



^a(Q high', Hi(Q high), H 2 (Q high)) = ^a(Q high] v l, v 2) 



(7.9) 



where v±, v 2 are the vevs of the EW minimum.^ In other words, to find the LHS of 
Eq. (|7.9p we simply integrate the RGE from Mx to this new Qhigh using as initial conditions 



In general the value of a running parameter depend on the current field point due to the field 
dependent thresholds involved in RGEs. 
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at Mx the outcome of the iteration procedure described in Sec. V. Evolving this set of 
values {Hi(Q h i gh ), H 2 (Q high), ^a(Q high, v 1 ,v 2 )} to Q*, using field dependent thresholds, the 
effective potential at the current field point can be constructed. 

VIII. CONCLUSIONS 

In the framework of MSSM, the fact that the top quark is heavy suggests an interesting 
possibility for explaining the spontaneous symmetry breaking at the EW scale, i.e., the 
radiative breaking senario. The key method to analyze such a senario is based on the RG 
equation. In describing the radiative symmetry breaking, the most primitive approach is 
to use the tree level Higgs potential with the RG running masses and couplings inserted. 
There exists, however, a serious technical problem in finding symmetry breaking solutions. 
Namely, the results often depend badly on the choise of the renormalization point Q at 
which the RG running is terminated, a fact that clearly reduces their reliability. 



As emphasized by the authors of Refs. |20|j2^ - p0| , |37| , addition of the one- loop corrections 



to the classical potential ameliorates the scale dependence at least for low energies. 
However, we have ascertained with surprise, that a careless treatment of these corrections 
for high energies formally jeopardizes the stability of the theory. Along a special direction 
(H-diag), where the magnitude of tree level potential is small, these corrections when 
carelessly treated predominate and produce an undesirable UFB escape. The reason behind 
such a failure must be sought in the inadequacy of a mass independent renormalization 
scheme (DR) to treat the very many mass scales of MSSM as one.0 The problem is that 
the decoupling of the various particles is not automatically included in the formalism and 
has to be incorporated. 

So to achieve our purpose we have tried to implement the decoupling theorem in a 
manner proposed by the authors of Ref. A simple way to incorporate it in the MSSM 



case is to treat the various particle thresholds as steps in the /^-functions, as well as in the 
one loop corrections of the scalar potential. Each time we cross a threshold the /3-function 
changes indicating that we have switched to a new effective field theory with the heaviest 
particle integrated out. At the same time, the associated particle's contribution to is 
dropped realizing in this way the process of decoupling. We stress here the role played by 
the renormalization scale choice, as given in Sec. VII. It should be wisely chosen in order to 
eliminate heavy particles whose participation induces a fake instability, while at the same 
time improve the convergence of the loop expansion (i.e., |V^| IV^ -*!). 

Employing the framework just stated and using the "merlin" minimization program |38| 
(for details see Appendix D) we have scanned the dangerous H-diag direction in its entirety 
in order to reveal unexpected local minima different than the true. This procedure has been 
repeated for a representative set of initial conditions at Mx, but the outcome was always 
the same: a bounded from below potential with a SU(3)c x U(l)y symmetric minimum at 
EW scale. 



20 For the case at hand (MSSM), the employnment of a mass dependent renormalization scheme 
would be an arduous task (for an application to the simple Yukawa model see ||). 
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Clearly, since the described method allows to undertake faraway excursions in the field 
space it can be utilized in less investigated situations. For instance, in the past various au- 
thors have focused on conditions (involving soft trilinear scalar couplings), which are needed 
to ensure that a particular SUSY model avoids an electric and/or color charge breaking 
ground state. However, to our opinion a careful numerical analysis of the impact the one- 
loop corrections would have on these matters, is still required. Due to their inherent 
complexity and for the sake of presenting analytic expressions, one usually resorts to getting 
round contrivances in order to efficiently deal with the problem. On the other hand, from our 
numerical point of view, we can directly attack the problem of radiative corrections at the 
cost of loosing track of analytic elegance. These issues will be the subject of a forthcoming 
publication |39 . 
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APPENDIX A: FIELD DEPENDENT SPECTRUM OF MSSM 



We give here all the field dependent mass eigenstates of the MSSM in the presence of a 
neutral Higgs background and the necessary expressions to compute £* of Eq. ([4.31). 



Similar formulae can also be found in Ref. f4T) |. To simplify somehow the expressions we 



switch to the notation Hi = x 1; 7i 2 = x 2 . The function f(x), frequently used below, has 
already been defined in Sec. II, as well as the rest of the notation. 

• Gauge Bosons (W ± ,Z) 



M\W±) = f (*? + xl) M\Z) = 9 —^(xl + x\) (A.1) 



and 



MW^ = MW±) = ^(M 2 ^)) (A.2) 

Mz) = Mz) = 3{92 J 2 9l) f(M 2 (z)] 



Leptons (£ : e, /i, r) 



(A.3) 
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and 



Y} 



8tt 2 



* 2 (^) = 



(A.4) 



Up-Quarks (Q:u,c,t) 



M 2 (Q) = (Y Q x 2 f 



(A.5) 



and 



= 



MQ) = -f* /(m 2 (Q)) 



(A.6) 



Down-Quarks (g : d, s, 6) 



M 2 ( (Z ) = (y s x 1 ) s 



(A.7) 



and 



oy2 



* 2 (?) = o 



(A.8) 



Gluinos 



M 2 (G) = Mf 



=* 2 (G) = 



(A.9) 



However, with a "vacuum energy" subtraction their contribution is canceled from 

Vl-loop- 



Charginos (Ci, C 2 ) 

Field dependent mass matrix: 



M| + -£ 2 (/^i + M 2 x 2 ) \ 

-g 2 {\xx x + M 2 x 2 ) n 2 + glx\ ) 



(A.10) 



Its eigenvalues M 2 (Ck) = \k(xi,x 2 ), k — 1, 2 are roots of the characteristic equation : 

\ 2 -b 1 \ + b = (A.ll) 

with hi = fi 2 + Mf + g\{x\ + x\) and 6 = (/^M 2 - g%x 1 x 2 ) 2 - Obviously b , bi > 
so there are two real and positive roots. Easily : 



cc ^ dec ^ 



+ 



61 <%i _ 2_dbo 



Xi dxi ' sign(A fe - A^ fe )P c 



(i,^= 1,2) 



(A.12) 
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where 



V c = yjb 2 - 46 = + M 2 y + g 2 {x x - x 2 f - M 2 ) 2 + g 2 [x x + x 2 f (A.13) 



Finally 



*i(Cfc) 



1 1 dA fc 
167r 2 a;,- 9a;,- 



/(A* 



(A.14) 



Neutralinos (JV 1? JV 2 , N 3 , JV 4 ) 
Field dependent mass matrix: 



M 



N 



( M 1 gxjV2 -gx 2 /V2\ 

M 2 -g 2 xi/V2 g 2 x 2 /V2 

gxi/\/2 -g 2 xi/\/2 -\i 

\-gx 2 /V2 g 2 x 2 /V2 -y, / 



(A.15) 



Mass eigenvalues M(Nk) = \k(xi,x 2 ), k = 1,2,3,4 are roots of the characteristic 
equation : 



A 4 + A 1 X 3 + A 2 \ 2 + A 3 X + A A = 



(A.16) 



It is straightforward to show that 
A x = -Mi - M 2 

A 2 = -/i 2 + M 1 M 2 - (g 2 + gl){x\ + x 2 2 )/2 

A 3 = /i 2 (Mx + M 2 ) - n{g 2 + gl)x 1 x 2 + (M ig 2 + M 2 g 2 )(x\ + x\)/2 
A A = -y 2 M x M 2 + n(M ig 2 + M 2 g 2 ) Xl x 2 



(A.17) 



Differentiation of Eq. ([A.16Q gives 



oc 2 doc ^ 



tk 



1,2) 



1 dA, 
— 7) — k 

OC j CJOC q 



ldA 2 ldA 3 ldA 4 
A fc H ^ — H — A fc + 



oc '} doc j 



OC j QOCrj 



oc j doc j 



A\\ + 3A.X 2 + 2A 2 \ k + A 3 



(A.18) 



Finally 



1 1 d\ 



327r 2 x,- dxi 



(A.19) 
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Sneutrinos (vg) (£ : e, //, r) 



M 2 (z>,) = m| + X - (g 2 + gl){x\ - x\) 



(A.20) 



and 



*i(£^) = -# 2 (^) 



ff 2 + ffl 
64tt 2 



(A.21) 



Other scalars (tpi,tp2) All other scalars can be grouped into independent subsets. 
These are: Sleptons (£i,£ 2 ), Up-Squarks (Qi, <2 2 ), Down-Squarks {qi,q 2 ) and Higgses 
(H+,H~), (h,H), (0 l5 2 ). 

The field dependent mass matrix, for each one subset, has the generic form 



m: 



P Q 
Q R 



(A.22) 



where tp runs on the previously mentioned pairs and the necessary quantities for each 
scalar subset are given in Table II and III. In order to find its mass eigenvalues 
M 2 ({pk) = \k(xi,x 2 ), k = 1,2 we have to solve the characteristic equation 



A 2 - (P + R)\ + PR-Q 2 







(A.23) 



Differentiation of the above equation gives 
1 d\ k 



i o{r + H) i x . dx . x . dx . 



X i 



dxi 



+ 



sign(A fe - Xt^yDif 



M = l,2) 



(A.24) 



where 



(P - R) 2 + 4Q 2 



(A.25) 



while the radiative corrections become 

1 



64tt 2 



(A.26) 



APPENDIX B: MASS EIGENSTATES AT SPECIFIC DIRECTIONS 

We present below mass eigenvalues and relevant quantities mentioned in Sec. III. Note 
that Yukawa couplings of the 3 rd family are denoted, as usual, by (Y t) Y b , Y T ) = (F u 3 , Y" b 3 , F e 3 ) 
and M s is the characteristic scale of SUSY breaking, which is typically an average of the 
relevant soft parameters. (Consult also Sec. II and Appendix A). 
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1. Direction (Hi,H 2 = Hi) = (r >, 9 = it/4) 
Gauge Bosons (5 = 1) 



2 2 2 2 

M 2( Z) = 9_±92 r 2 M 2 [W ± ) = 9_2]_ (R1) 



Fermions (S = 1/2) 



M 2 (t) = - Y t 2 r 2 M 2 (b) = ~ >fr 2 M 2 (r) = 1 y-', 2 

2 2 

M 2 (Ci) = ^f- M 2 (C 2 ) 



1 




2 










c/2 r2 




2 



(B.2) 



For the neutralino eigenproblem we have assumed a perturbative solution [41[ of the 
form 

M(Ni) = \i (i = 1,2,3,4) with 

\l = A 2 = 6 r ^3 = & A 4 = £4 

Since Aj are roots of the characteristic equation A 4 + A 1 X 3 + A 2 X 2 + A 3 X + A A = (see 
Appendix A), we know from elementary algebra that A, satisfy the following system 
of equations 

Ai + A 2 + A 3 + A 4 = -A 1 
A1A2 + A2A3 + A3A4 + A4A1 + A1A3 + A2A4 = A2 
A1A2A3 + A2A3A4 + AiA 2 A4 + A1A3A4 = —A3 

AiA 2 AsA4 = A4 



By comparing the highest order terms in r we imediately conclude that 

t e 1 92 ± 92 

whereas £ 3 , £ 4 are solutions of 

-g 2 e + {x- w 2 )i + = 0, 

g 2 =g 2 + g 2 and X = M l9 2 + M 2 g 2 
So the relevant eigenvalues become 



M 2 (N 1 ) = M 2 (N 2 ) = g2 + 92 r 2 M 2 (N 3 )=0(M 2 ) M 2 (N A ) = 0(M 2 S ) (B.3) 
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• Higgses (S = 0) Solving the eigenvalue problem to highest order in r we get 
M 2 (H + ) = ^— M 2 (H-) 



M 2 (H) = U ' ™ r 2 M 2 (h) 



2 2g 2 r* 
9 2 + 9 2 2 w2/ fc x _ m l + m l ( m l ~ m l) 2 (B.4) 



2 ' ~ 2 2(^ + ^ )r2 

M 2 (0!) = 0(M 2 ) M 2 (0 2 ) = C»(M|) 



Super scalars (S = 0) 



M 2 (z> 3 ) = 


0(M 2 S ) 


M 2 {f 1 ) = 


\yy 


M 2 {h) = 


\yy 

2 * 


M 2 (t 2 ) = 


\yy 


M 2 (bi) = 


2 b 


M 2 (b 2 ) = 


\yy 



M 2 (f 2 ) = 1 -Yy 



2. Direction (0,W 2 ) = (r », = vr/2) 
Gauge Bosons (S* = 1) 



(B.5) 



2 2 2 2 



• Fermions (S = 1/2) 

M 2 {t) = Y 2 r 2 M 2 {b) = M 2 (r) = 

/y 2 M 2 ( R7 ) 

92 r 

For the neutralino sector as before, we assume a perturbative solution now with a 
slightly different form 

Ai = £ir A 2 = 6r A 3 = £ A 4 = £/r 2 
Solving to highest order we obtain 



c _ c \ g 2 + 9 2 2 c , _ Migj + M 2 g 2 2n 2 M 1 M 2 

~K\ - ?2 - V o ^3 - o T o ^4 



2 ^ + ^ ~ M l9 l + M 2 g* 

and the relevant eigenvalues become 

M 2 (iV 1 ) = M 2 (iV 2 ) = g2 + 92 r 2 M 2 (N 3 ) = 0(M 2 ) M 2 (N 4 ) = % (B.8) 

2 v 
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Higgses (S = 0) 

M\H+) = g —±A T 2 M\H-) = Z9l±A r 2 

M 2 {H) = 3{9 2 + 9l) r 2 M 2 {h) = -t±& r 2 (B.9) 



M 2 ((pi) = -t±A r i m 2 (0 2 ) = g ^±A r 



Super scalars (S = 0) 

M 2 (z> 3 ) = 



9 2 


+ 92 




4 






g 2 r 2 




2 








yy 


9 2 




3 






g 2 r 2 




6 





.U' J (7,)= ^— M 2 (f 2 )= g!±M r 2 



..uV, i = v;-/- 2 ^ M 2 (t 2 ) = K t V + 9 * 3 ^ 2 r 2 



(B.10) 



12 



^(60 = ^- M 2 (6 2 ) = ^^r 2 

Finally, the expressions (3.2) for the one-loop effective potential along (0, W 2 ) direction 
require the following quantities 

TT - 9 2 +9i rr _ 92 rr _ 2 „_ ^ + ^1 

Ui - — 2 — T 2 — 4 

£ 2 + 3s 2 2 rr _3((? 2 + ^ 2 2 ) rr # 2 rr _ v2 Q 2 -3g 2 



(B.12) 



8 12 9 J ~ iU ~ r 3 ~ iX ~ r ' 12 



and 



u 3 = 


9l 




u 4 


u 7 = 


9 2 
6 






u w = 




3 




2, 4, 


6, 6, 


1, 6, 


6) 
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APPENDIX C: REALIZATION FOR /5-FUNCTION THRESHOLDS 

The generic form of a step function appearing in MSSM RGEs with one-loop thresholds 
is 

e . = e(Q 2 - Ti{Hi, H 2 )) (i = 1, 2, ... , 22) (C.l) 

where Tj defines the decoupling point. Following the notation of Sec. II and Appendix A, 
we shall present here an implementation of the various Tj we have used in our investigation. 
• Field configurations near EW breaking minimum (< 300 GeV) 
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Ti = M 2 


T 2 = 


Ml 


n-~ 




T 4 = m 2 Hl 




T 5 = ™h 2 


T 6 = /i 2 


T 7 = 




T 8 ~- 


2 

= m% 

W3 


T 9 = m | 2 




T ™ = m |x 


Tu = m} 

U 3 


Ti 2 = 


: m} 


Tl 3 


= ml { 


T u = nA c 

U 3 




T15 = m\ c 


Tw = m| c 


Tl 7 = 


9 

= mf 


Tis 


= mf 


T 19 = m 2 Li 




^20 = m 2 e, c 

a 3 


T 2 i = m\ c 


T 22 = 


= m| f 












Otherwise?]] 
















Ti = \M 2 (N 3 )\ 


T 2 


= |M 2 (C 2 ) 


1 


T 3 = 




T 4 ~- 


= \M 2 (<Pi)\ 


T 5 = \M 2 (H+)\ 




= |M 2 (iV 4 ) 


1 


T 7 = 


M 2 (Ci)\ 


Ts~- 


= \M 2 (ti)\ 


T 9 = \M\di)\ 


Tio 


= |M 2 (^ 


)| 


Tn = 


\M 2 (t 2 )\ 


Tn 


= \M 2 (5 2 )\ 


T n = \M 2 (u 2 )\ 


T u 


= \M 2 (bi) 


1 


Ti 5 = 


\M 2 (Si)\ 


Tie 


= \M 2 (di)\ 


Ti 7 =\M\fi)\ 


Tis 


= WXvi 


)l 


T w = 


\M 2 {ei)\ 


T 2 q 


= \M 2 (f 2 )\ 


T 2 i = \M 2 (fi 2 )\ 


T 22 


= \M\e 2 ) 


1 











(C.2) 



(C.3) 



APPENDIX D: NUMERICAL ISSUES 

• Rounding error: Squeezing infinitely many real numbers into a finite number of bits 
(binary digits) requires an approximate representation. Given any fixed number of bits, (e.g., 
Double Precision has 64 bits) most calculations with real numbers will produce quantities 
that can not be exactly represented using that many bits. Therefore, the result of a floating- 
point calculation must often be rounded in order to fit back into its finite representation. 
This rounding error is the characteristic feature of floating-point computation. 

Estimation of this kind of error is in general an extremely complicated process. However, 
for the trivial operations (+, — , *, /) one can give bounds for the rounding errors involved 
p2] , |43[| , From our perspective, the computation of in case of huge field dependent 

(FD) mass eigenvalues, might suffer from large rounding error. This error could arise from 
subtractions of numbers that are nearly equal or additions and subtractions of numbers that 
differ greatly in magnitude. Clearly, if it is greater in absolute value than V^°\ then our 



resulted Vi-i oop would be untrustworthy. One can show |43| that a crude upper bound for 
the relevant rounding error due to summation is 

E 8 (V^)<5^\vi\ 8 = ^ m 

For Double Precision accuracy t = 16 
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We use M 2 (V) to denote the field dependent mass eigenstates of the particle V . 
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Fortunately enough, one can reduce rounding error in these cases, by grouping the operands 
according to relative size, so that as much as possible operations will be performed between 
numbers with similar magnitudes. The recommended method is sorting the partial sums 
(V^) with respect to their absolute value before summing them.0 Afterwards, summation 
is performed by using the Kahan formula [ 4"4"| . Consequently, we have to restrict ourselves 



only to those field configurations for which 

\V^\ > E S (V {1) ) (D.2) 

is satisfied. Indeed, from our numerical study it seems that for all \Hi\, \H.2 \ < 10 11 GeV 
Eq. ( p.2|) holds. So one should not worry for this kind of error. 

• Catastrophic canselation error: The evaluation of any expression, containing an ad- 
dition of quantities with opposite signs, could result in a relative error so large that all the 
digits are meaningless. When subtracting nearby quantities the most significant digits in 
the operands match and cancel each other. 

Especially, a catastrophic canselation occurs when the operands are approximately equal 
and subject to rounding errors. However, a formula that exhibits catastrophic canselation 
can sometimes be arranged to eliminate the problem. We have faced this problem in com- 
puting some FD mass eigenstates of various MSSM sectors. These, for most MSSM scalar 
particles, are determined by solving the eigenvalue problem of a generic 2x2 mass matrix 

m *=(qp) < d - 3 > 

The solution of its characteristic equation 

x 2 + bx + c = 



b = P + Q, c = PQ-R 2 
involves the computation of two "dangerous" quantities 



(D.4) 



A = b 2 - 4c and Z± = -b± Vb 2 - 4c (D.5) 

When b 2 ~ 4c catastrophic canselationP] can cause many of the accurate digits to disap- 
pear leaving behind mainly digits contaminated by rounding error.0 Luckily enough, this 



This procedure gives the smallest upper bound in Eg, but it does not necessarily give the 
smallest error in practice. 

23 The quantities b 2 and 4c are subject to rounding errors since they are the results of floating 
point multiplications. 

24 Indeed we found that in some cases round off error completely ruins the results by turning A 
into an imaginary quantity. 
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shortcoming can be fixed by rearranging the terms of A. Namely, we can use the following 
expression 



A = (P + Qf - 4(PQ - R 2 ) = {P- Qf 



AR 



(D.6) 



everywhere in our computations. On the other hand, if c is small then one of Z± will involve 
the subtraction of b from a very nearly equal quantity (A) and the associated root comes 
out with large inaccuracy. However, a definitely more accurate way to compute the roots of 
( [Dip is found in [g§ 



xi = q 



b + sign (6) V^ 2 — 4c 
c 
9 



and 



(D.7) 



X-2 



The expression x 2 — y 2 is another formula that exhibits catastrophic cancelation (when 
x ~ y). To compute x 2 — y 2 , as accurate as possible, we apply |f42] 



(x + j/)(x-y), 
— 



if 75 < 
otherwise 



< V3 



(D.8) 



• Diagonalization: The optimum strategy for finding eigenvalues is first to reduce the 
matrix to a simple form and afterwards begin an iterative procedure. The most efficient 
program for finding all eigenvalues of a symmetric matrix, is a combination of the House- 
holder tridiagonalization and the QL algorithm with implicit shifts (for details see [f45| , [46|1 ). 
Note that when dealing with a matrix, whose elements vary over many orders of magnitude, 
it is important that the smallest elements are in the top left-hand corner. This is because 
the reduction is performed starting from the bottom right-hand corner and a mixture of 
small and large elements there can lead to considerable rounding errors. One possible way 
to overcome this difficulty is to perform a trivial reflection through secondary diagonal when 
matrix elements are not properly aligned. Of course this is an orthogonal transformation so 
the eigenvalues are not affected. 



• Minimization aspects: Minimization of the effective potential Vi-i oop was performed by 
running the "merlin" software package on a SGI Origin 2000 supercomputer. Merlin is an 
integrated environment designed to solve optimization problems of the following form : 

Find a local minimum of the function 

/(£), x e H N , x = [x 1 ,x 2 ,...,x N ] T 



under the conditions 

Xi e [ai,bi] for i = 1,2, ...,N 
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It contains implementations of powerful minimization algorithms. In particular, there are 
two direct methods (ROLL, SIMPLEX) that use no derivative information and three algo- 
rithms from the conjugate gradient family (Fletcher-Reevs, Polar-Ribiere and the generalized 
Polar- Ribiere) . Also from the Quasi- Newton methods the DFP and several versions of BFGS 
method are coded. 

To improve minimization's effectiveness, we have made a linear change of independent 
variables (scaling) so that the values of the new variables at the minimum are of order unity 



47| . For convenience, we have also imposed an upper bound on the domain of Vi„/ oop coding 



in our routines the following potential functionP] 

—loop J 

if \Hi\, |^2 1 < 5*00 m q n 

SGI ~ \ V Big otherwise (D ' 9) 

We have tried many minimization algorithms to our potential function. Of course, due 
to lack of analytic derivatives and the extremely "flat" behaviour of (p.9|) , near H-diag 
direction, direct methods have an advantage over the gradient ones. Indeed, by means of 
the SIMPLEX algorithm we were able to scan the H-diag direction in its entierty for probable 
local minima. 



Our choise is ~ 10 11 GeV and V Big = 10 300 (GeV) 4 . 
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TABLES 



TABLE I. Experimental bounds on super symmetric particles (in units of GeV). 



Neutralinos 


M Nl 


> 


10.9 


M N2 > 45.3 


M Ns 


> 75.8 




> 


127 


Charginos 








M Cl > 65.7 


M c 


2 > 99 








Sneutrino-Higgses 




> 


43.1 


M h > 62.5 


M H 


> 77.5 


M H± 


> 


54.5 


Squarks-Sleptons 




> 


224 


Mg > 58 




> 55.6 


Mr 


> 


45 



TABLE II. Degrees of freedom for all MSSM mass eigenstates in a neutral Higgs background. 



Eij 


$enstate 


Spin (S) 


Color (C) 


Helicity (M) 


Gauge Bosons 


(W±) 


1 


1 


2 




(Z) 


1 


1 


1 


Leptons 


(£ : e,n,r) 


1/2 


1 


2 


Up-Quarks 


{Q:u,c,t) 


1/2 


3 


2 


Down-Quarks 


(q : d,s,b) 


1/2 


3 


2 


Gluinos 


1/2 


8 


1 


Charginos 


(Ci,C 2 ) 


1/2 


1 


2 


Neutralinos 




1/2 


1 


1 


Sneutrinos 


in) 





1 


2 


Sleptons 


{hM 





1 


2 


Up- S quarks 


(Qi,Q 2 ) 





3 


2 


Down-Squarks 


(91,92) 





3 


2 


Higgses 


(H+,H~) 





1 


2 




(h,H) 





1 


1 




(01,02) 





1 


1 



TABLE III. Mass matrix entries for some MSSM eigenstates in a neutral Higgs background. 



Eigenstate 


P 


Q 


R 


Sleptons 


^|c + ^l + T(^2-^l) 


Ye(Axi +11x2) 


m) +Y?x\ + g2 ~ 9 Hxi x 2 2 ) 


Up-Squarks 

(Qi,Q 2 ) 




-Yg(fJ,Xi + Ax 2 ) 


m 2 >+Y 2 x 2 + 9^9i {x 2_ x 2 ) 


Down-Squarks 
(91,92) 


+ + £(*!-*?) 


Y q (Ax 1 + jix 2 ) 


ml + Y^xl + ^ixl-xl) 


Higgses 
(H+,H-) 
(h,H) 

(01,02) 


ml + g2+ / 2 {^l x\) 


-mg + -fxix 2 

2 a 2 +g 2 

-ml 


ml + 9 Y 2 {xl + 
m 2 2 + g2+ /i(3x 2 2 x\) 
m| + s > 2 (x| x\) 
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FIGURES 




FIG. 1. The 1-loop potential as a function of the "generalized radious" R g {4>i) = y</>i + 4>\ 
shown in the A = -400 GeV, M = 60 GeV, M 1/2 = 200 GeV and tan/3 = 2 case. {<f> u fcifa)) 
the lowest point of H-Diag direction at an arbitrary field value 4>\ = Ti.\{Qhigh)- 
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FIG. 2a. Decoupling process of the top-stops field dependent mass eigenvalues as we are 
moving along the line x\ + X2 = (pi + ^2(^1)5 where Xi = TLi. ((pi, <fe(</>i)) represents the lowest 
point of H-Diag direction when <p\ = Tii(Qhi g h) — 2 x 10 3 TeV (initial conditions at Mx same as 
in Fig. 1). 
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FIG. 2b. Same as Fig. 2a but now for charginos and bottom-sbottoms. 




FIG. 3a. Plot of the u(x) used in the definition of the critical scale Q*. 




FIG. 3b. The renormalization scale used for RG improvement of the effective potential as a 
function of x = \ log . 



